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^ ■ Context. The 2008 outburst of the atoll source IGR J17473-2721 was observed by INTEGRAL, RXTE and Swift. Tens 

of type-I X-ray bursts were found in this outburst. 

Aims. Joint observations provide sufficient data to look into the behavior of IGR J17473-2721 at the rising part of the 
outburst. The relation between the duration of the bursts and the accretion rate and the nature of the corona producing 
the observed power-law component can therefore be studied in detail. 
, ^ Methods. We analyze observational data of IGR J17473-2721, focusing on the spectral evolution during the state 

, transition from quiescent to low hard state (LHS), and on the flux dependence of the type-I X-ray bursts along the 

outburst. 

Results. We find that the joint energy spectra can be well fitted with a model composed of a blackbody and a cutoff 
power-law, with a cutoff energy decreasing from ~ 150 keV to ~ 40 keV as the source leaves the quiescent state toward 
^ , the low hard state. This fits into a scenario in which the corona is cooled by the soft X-rays along the outburst evolution, 

^ ■ as observed in several other atoll sources. By using the flux measured in the 1.5-30 keV band of the type-I bursts during 

' the outburst, we find that the linear relationship between the burst duration and the flux still holds for those bursts 

I that occur at the decaying part of the low hard state, but with a different slope than the overall one that was estimated 

with the bursts happening in the whole extent of, and for the rest of the low hard state. The signiflcance of such a 
dichotomy in the type-I X-ray bursts is ~ 3 a under an F-test. Similar results are hinted at as well with the broader 
energy-band that was adopted recently. This dichotomy may be understood in a scenario where part of the accreting 
^ ' material forms a corona on the way of falling onto the surface of the neutron star during the decaying part of the low 

00 , hard state. Based on the accretion rates of the preceding LHS, estimated from type-I X-ray bursts and from persistent 

emission, at least for IGR J17473-2721, most of the accretion material may fall on the neutron star (NS) surface in the 
' LHS. Considering the burst behavior in the context of the outburst indicates a corona formed on top of the disk rather 

' than on the NS surface. 

Key words, star: neutron - individual: IGR J17473-2721 -X-rays: bursts 

I : 

• • . 1. Introduction In general, the spectra of LMXB (both BH and NS bina- 

' ^. 1 ■ ■ rTT,fT^r-r>\ , • • , ries) are described by two-component models, including a 

Low-mass X-ray binaries (LMXB) containing neutron stars ^oft/thermal (e.g., a blackbody) and a hard/Comptonized 

(NSs) are classified as atoll or Z sources according to then ^^^^^^^^^ ^ power-law shape with a cutoff at tens to 

^ different spectral and timing evolution during outburst j^^^^reds of keV). It is believed that the hard component 

■ (Hasmger and van der Klis 1989). Atoll sources have a typ- , • +i, • ri j. i-i- ■ t ti-uu i 

. . n 1 1- 1 evolves m the inverse Uompton scattering of sott/thermal 

ical luminosity < 10% Eddmgton, and mostly evolve ^^^^^^^ electrons in the corona. With respect 

from the island to the banana state, analogous to that of ^^^^^^ -^^^^^^ ^^-^ uncertainties remain for 

the black hole (BH) counterparts which evolve from the ^^^^ ^^^^^ ^^^^ ^^^^^ electrons in 

LHS to the high/soft states (HSS). Z sources are usually ^^^^^^ energized up to ~ 100 keV, and what is the 

brighter, with softer spectra and weaker variability, m anal- , j i j-- /'r> u p m i onno\ 

f +i, w«Q • "nw • corona geometry and location (Reynolds & Nowak 2003). 

ogy to the HSS m BH b inaries. Regarding the latter, the Eastern (Mitsuda et al. 1989) and 

Send offprint requests to: Yu-Peng Chen the Western models (White et al. 1988) have been proposed 

* Based on observations with INTEGRAL, an ESA project ^or NS binaries; whereas models with the so-called "sand- 

with instruments and science data centre funded by ESA wich" geometry, "sphere-^disk geometries", and "patchy 

member states (especially the PI countries: Denmark, France, corona" were proposed as well for BH binaries (see Figure 6 

Germany, Italy, Switzerland, Spain), Poland and with the par- of Reynolds & Nowak 2003). The Eastern model considers 

ticipation of Russia and the USA. that the thermal emission arises from the disk and that the 
Correspondence to: chenyp@mail.ihep.ac.cn 
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Comptonization arises from hot plasma around the neutron 
star. Instead, the Western model considers that the ther- 
mal emission arises from the boundary layer and that the 
Comptonization arises from the disk. 

For the Comptonized components of both BH and NS 
binaries in LHS, the high-energy cutoff is present at tens to 
hundreds of keV. This can be an indicator of the electron 
Maxwellian energy distribution in the corona and approx- 
imately twice the electron temperature (Hua & Titarchuk 
1995). Studying the evolution of the cutoff energy during 
outbursts is essential for understanding the accretion pro- 
cesses. For the LHS of BH binaries, the power-law is cut at 
~ 100 keV, while no cutoff was detected up to 10 MeV in 
the HSS of seven transient BH candidates: GRO J0422-h32, 
GX 339-4, GRS 1716-249, GRS 1009-45, 4U 1543-47, 
GRO J1655-40, and GRS 1915+105 (Grove et al. 1998). 
Recently, for the well-studied BH binary GX 339-4, the cut- 
off energy monotonically decreased, starting from '^200 keV 
and reaching ~60 keV, while the flux of the whole X-ray 
band increased during the LHS before the HSS (Motta et 
al. 2009). Interestingly, the cutoff energy is still detected 
(at ^ 100 keV) in high-intermediate and soft-intermediate 
states and suddenly disappears in the HSS state. For the NS 
binaries, no cutoffs were observed in the bright Z sources 
(see the review by Done et al. 2007). For some of the atoll 
sources observed, the cutoff energy is seen at more than 20 
keV in LHS (Barret & Vedrenne 1994). For instance, Ghen 
et al. (2006) found that the cutoff energy changed from ~ 
3 keV to > 30 keV during the spectral transition from the 
HSS to the LHS of 4U 1608-522, proposing that the exis- 
tence of a solid surface or the influence of a strong magnetic 
field may have a great impact on this difference with BH 
binaries. 

However, the changes in the cutoff energy of NS bina- 
ries while the source transitions from quiescence to the LHS 
are much less addressed, partially because the fact that the 
outbursts are usually caught when the source leaves the qui- 
escent state where the emission dominated by hard X-rays 
is too weak. The joint observations of the initial rising phase 
of the outburst from IGR J17473-2721 by INTEGRAL and 
RXTE provide a nice opportunity to conduct this study. 
The atoll source IGR J17473-2721 has shown a variety 
of spectral states/transitions in the 2008 outburst. At the 
rising phase of this outburst, the source showed hard X- 
ray flares and stayed in the low/hard state for two months 
prior to transition to the high/soft state, forming the so- 
called hysteresis. IGR J17473-2721 is thus the third NS 
binary with hysteresis after 4U 1908-1-005 (Aql X-1) and 
4U 1608-522 (e.g. Gladstone et al. 2007). The outburst with 
a low/hard state prior to a high/soft state resemble Aql X-1 
and XTE J1550-564. The long-lived preceding LHS makes 
IGR J17473-2721 to resemble the behavior of outbursts 
seen in black hole X-ray binaries like GX 339-4. Zhang et 
al. (2009) presented details of RXTE observations on the 
outburst, but data for the initial part were not yet available. 

Type-I X-ray bursts manifest themselves as a sudden 
increase (typically by a factor of 10 or greater) in the X- 
ray luminosity of neutron star systems (NSs), with ~90 X- 
ray bursters detected in our Galaxy (Grindlay et al. 1976, 
Behan et al. 1976). Type-I bursts are caused by unstable 
burning of accreted H/He on the surface of neutron stars 
in low-mass X-ray binary (LMXB) systems, in contrast to 
type-II bursts, that are thought to be caused by accretion 
instability (for reviews, see Lewin et al. 1993, Gumming 
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Fig. 1. Significance map of IGR J17473-2721 during the 
outburst in 2008, as observed by INTEGRAL. The contour 
start at a significance level of 5cr with steps of 50cr. 

2004, Strohmayer & Bildsten 2006, and Galloway et al. 
2008). The properties of type-I X-ray bursts depend on 
the composition of the burning material, how stable the 
burning of that fuel is, and if there is any left-over fuel 
from previous bursts. Several other factors can also signifi- 
cantly change the burst properties: individual element reac- 
tion rates, fuel mixing in the burning layer, inhomogeneous 
distribution of fuel, and other aspects. 

Tens of type-I bursts were detected in IGR J 174 73- 
2721 during outburst by RXTE, INTEGRAL, Swift, and 
SuperAGILE (Altamirano et al. 2008, Galloway et al. 2008, 
Ghen et al. 2010, Ghenevez et al. 2010). Chen et al. (2010), 
who found 16 bursts, estimated that the distance is likely 
to be 6.4±0.9 kpc based on the peak flux of three bursts 
that show photospheric radius expansion. The duration of 
the 16 type-I bursts occurred in the 2008 outburst is found 
to correlate with the Eddington ratio and with two paral- 
lel evolution groups (Chen et al. 2010). Analysis of more 
complete RTXE data by Ghenevez et al. (2010) revealed 
that 57 bursts (42 bursts among those detected by RXTE) 
were born in the 2008 outburst of IGR J17473-2721, with 
an overall correlation between the burst duration and the 
flux. However, Ghenevez et al. (2010) adopted an expanded 
energy band (0.1-200 keV) to calculate the bolometric flux 
of the persistent emission, far beyond RXTE^s best cov- 
erage (1.5-30 keV), which may result in flux uncertain- 
ties as large as ~40%. This stimulates us to re-investigate 
the duration/flux relation with this burst sample. Finally, 
INTEGRAL data on the initial rising part of the outburst 
were never reported before, which could provide clues on 
the formation of the outburst from quiescence and help in- 
vestigating the initial cooling process of the corona. We 
report on this here as well. 

2. Observations and data analysis 

The available INTEGRAL observations, when IGR J17473- 
2721 had an offset angle less than 14° during the outburst 
in 2008, comprised about 463 SCWs, adding up to an expo- 
sure time of 1.3 Ms (covering revolutions 0663-0723 -one 
revolution covers '^three days- between MJD 54545.1 and 
MJD 54750.0). However, data only scatter over the rising 
and decay phases in this outburst, and only the former is 
used here because of the low signiflcance of detection in the 
decay phase. The data reduction was performed using the 
OSA version 9.0. ISGRI and JEMX data are used to extract 
the spectra and lightcurves. All sources within the FoV that 
are brighter than, or comparable to, IGR J17473-2721 were 
taken into account in extracting the source spectrum. 
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Fig. 2. Upper panel: PCA (2-12 keV) lightcurve of IGR J17473-2721 with a time bin of 16 seconds, in which the spikes 
are type-I X-ray bursts. The color bar above the PCA lightcurve indicates the burst groups for the subsequent analysis 
of the relationship between the burst duration and the corresponding persistent flux. Middle panel: the BAT (15-50 keV) 
lightcurve of IGR J17473-2721 (black) and GRO J1750-27 (red) during the outburst in 2008. Lower panel: the ISGRI 
(20-200 keV) lightcurves of both sources monitored by INTEGRAL. In both of the lower panels the bin size is one day. 



Public data from RXTE (Gruber et al. 1996) on 
IGR J17473-2721 cover the entire 2008 outburst between 
April 2008 and October 2008, and include 182 RXTE /VGA 
pointed observations, with the identifier (OBSID) of pro- 
posal number (PN) 93064, 93093, and 93442 in the High 
Energy Astrophysics Science Archive Research Center 
(HEASARC). These observations add up to 475 ks of 
exposure time on the source, and scatter over the entire 
outburst. The analysis of the PCA data was performed by 
using HEAsoft v. 6.6. We filtered the data using the stan- 
dard RXTE /PCA criteria. Only the PCU2 (in the 0-4 num- 
bering scheme) was used for the analysis, because this PCU 
was 100% working on during the entire observations. The 
background file used in the analysis of PCA data is the most 
recent one for bright-sources found at the HEASARC web- 
sitcQ; the detector breakdowns have been removed. Data 
from cluster 1 of the HEXTE system were used to produce 
the higher energy spectra. 

An additional 1% systematic error is added to the spec- 
tra because of calibration uncertainties, if not otherwise 
specified. The X-ray spectra are fitted with XSPEC vl2.5.0 
and the model parameters are estimated with a 90% confi- 
dence level. In fitting the joint RXTE k INTEGRAL data 
points, the JEMX, HEXTE, and ISGRI spectra were nor- 
malized to the PCA one. 

There is only a single OBSID 00031182001 from Swift on 
IGR J 17473-2721 during the 2008 outburst, on 2008 March 
3 (MJD 54556). During the observation, the XRT data were 



^ pca_bkgd_cmbrightvle_eMv20051128.mdl 



performed with Windows Timing (WT) mode and Photon 
Counting (PC) mode. We extracted PC mode events using 
the tool XSELECT V2.4a as a part of HEASOFT v6.9 
package, and pile-up effects were corrected. 



3. Results 

3.1. Outburst 

We find that IGR J 17473-2721 has a detection significance 
of ■~182cr in the 20-200 keV band during the period between 
MJD 54545.1 and 54750.0 (see the INTEGRAL map in Fig. 
ID). The nearby source GRO J1750-27, that is 48.92 arcmin 
away from IGR J17473-2721, was out of the RXTE /PCA 
field of view during the observations. The field of view of 
RXTE/PCA is -1 degree, RXTE/PCA can't distinguish 
two sources separated by ~ 1 degree if one is located at 
the center of the field. RXTE/PCA can avoid contami- 
nation from GRO J1750-27 by pointing ^0.2 degree away 
from IGR J17473-2721. Fig.Ilshows the BAT0 and ISGRI 
lightcurves of both sources in bins of one day. 

We extracted the RXTE/PCA lightcurves in 2-12 keV 
with PCA standard 2 model and binned it into observation 
identifiers (OBSID), with the background removed. This 
is shown in the top panel of Fig. [31 where for comparison 
the INTEGRAL/lSGRl lightcurves in 20-200 keV are also 



^ See the Swift/BAT transient moni- 

tor results provided by the Swift Team at 
http://swift.gsfc.nasa.gov/docs/swift/results/transientsl 
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Unfolded Spectrum 



Fig. 3. Top two panels: RXTE /PC A lightcurves of IGR 
J17473-2721 (2-12 keV, time bin is fixed to one obsid) 
and INTEGRAL/lSGRl, 20-200 keV, time bin is fixed to 
one day) during the rising phase of its outburst in 2008. 
The temperature (T) and radius (R) of the blackbody, the 
photon index (F), and the cutoff energy resulting from fit 
results with Snnft /XRT (0.5-10 keV), RXTE /PC A (3-30 
keV), RXTE/EEXTE (30-100 keV), INTEGRAL/ JEMX 
(5-25 keV), and INTEGRAL /ISGRl (20-200 keV) are also 
given in the remaining four panels (the time bin of the bot- 
tom panels is one revolution, roughly happening every three 
days). 



shown. From Fig. [3] it is possible to see that IGR J17473- 
2721 left from quiescent state and started to be active at 
around MJD 54555, with the flux increasing in both en- 
ergy bands. Roughly 20 days later, at MJD ~54575, the 
source stayed at the LHS and the corresponding results are 
reported in Zhang et al. (2009). 

To analyze the spectral evolution of the initial rising 
phase of the outburst from IGR J17473-2721, the ISGRI 
data were divided into ten groups based on the revolu- 
tion number: rev. 663, 664, 667-674. For each group we 
fitted the spectra in the 20-200 keV band with a power- 
law, a cutoff power-law, or cutoff power-law, both with 
1% systematic error included. No RXTE or Swift ob- 
servations are available during the first two revolutions 
(663 and 664), and the source is too weak to have a 
meaningful JEMX spectrum. There is one Swift obser- 
vation quasi-simultaneous to the INTEGRAL revolution 
667, therefore the joint Swift /XRT, INTEGRAL/ JEMX 
and INTEGRAL /ISGRl was derived. IGR J17473-2721 
is outside the FoV of JEMX during rev. 668, the 
joint spectrum was obtained from INTEGRAL/lSGRl, 
RXTE /PGA and RXTE /REXTE observations. For rev. 




Fig. 4. Joint XRT (purple), JEMX (blue), ISGRI (green) 
spectra of rev. 667 of IGR J17473-2721 during the outburst 
in 2008, fitted with model: wabs*cutoff. 



669-674, we fitted the combined INTEGRAL/ JEMX, 
INTEGRAL/lSGRl, RXTE /PGA and RXTE /REXTE 
spectrum. RXTE OBSIDs where bursts occurred were ex- 
cluded from the spectral analysis. 

We extracted for each data group the spectra from the 
combined XRT data in 0.5-10 keV data, the PCU 2 data in 
the 3-30 keV band, the JEMX data in the 5-25 keV band, 
the HEXTE cluster 1 data in the 30-100 keV band and 
ISGRI data in the 20-200 keV band. We used the blackbody 
model (bbodyrad in Xspec) to fit the thermal component 
-because the disk emission is hardly visible owing to the 
relatively steep disk inclination angle (Zhang et al. 2009)- 
and a cutoff power law for the hard component fTable 13 . ip . 
The iron line is fixed at 6.4 keV and the absorption is set to 
3.8x10^^ atoms/cm^. A constant is introduced to account 
for the normalization difference among XRT, PGA, JEMX, 
HEXTE and ISGRI. 

We find that at the beginning of the outburst (rev. 663 
and 664 between MJD 54545 and 54550), and where only 
ISGRI spectra are available, a power law with a photon in- 
dex of ~ 2.2 (typical error of 0.2-1.0) can fit the spectra 
in the 20-200 keV band well, with reduced 1-6, 1.2 re- 
spectively (d.o.f. is 7, 7). For rev. 667 (MJD 54556), the 
joint XRT/ISGRI/JEMX data can be fitted with an ab- 
sorbed cutoff power law model, resulting in a Xred (l-H 
for 134 d.o.f.) comparable to that derived with a simple 
power law model (1.15 for 135 d.o.f.) (Fig.H Table[Sl]). By 
adding a blackbody and a Gaussian line fixed at 6.4 keV 
as additional components, the fit can be slightly improved 
(Xred/d-o-f- ^1.04/130), but the F-test shows a probability 
of only 0.022. 

From rev. 668 (MJD 54559) onward, the thermal com- 
ponent can be satisfactorily fitted with a blackbody model, 
while the fit of the hard component by a simple power law 
turns worse, with a reduced larger than 2 fTable [5TT|) . 
Adopting the hard component cutoff power law re- 

sults in better fits, with reduced improving to 1 (e.g.. 
Fig. [5] shows the spectral fit during rev. 674 as an exam- 
ple) . During the initial rising phase of the outburst we find 
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Unfolded Spectrum 




Fig. 5. Joint PCA (black), JEMX (blue), HEXTE (red), 
ISGRI (green) spectra of rev. 674 (the last revolution 
of INTEGRAL in the rising phase) of IGR J17473- 
2721 during the outburst in 2008, fitted with model 
wabs* (bbodyrad+cutoff+Gauss) . 

Table 1. The spectral fit results of the outburst with dif- 
ferent models. 



rev 


powerlaw 


cutoff 


cutoff+bbodyrad+gauss 






Xrcc 


(d.o.f) 


663" 


1.6(7) 


1.96(6) 




664" 


1.2(7) 


1.45(6) 




667" 


1.15(135) 


1.11(134) 


1.04(130) 


668'= 


5.22(97) 


3.91(96) 


1.30(92) 


669" 


7.33(134) 


3.45(133) 


1.04(129) 


670'' 


4.67(138) 


2.95(137) 


0.67(133) 


671'' 


6.26(138) 


2.39(137) 


0.68(133) 


672'' 


10.1(137) 


4.10(136) 


0.84(132) 


673'' 


12.9(138) 


3.97(137) 


0.71(133) 


674'' 


11.1(135) 


3.22(134) 


0.83(130) 



" with ISGRI data 

" with ISGRI+JEMX+XRT data 

with ISGRI+PCA+HEXTE data 
'' with ISGRI+JEMX+PCA+HEXTE data 



that the temperature of the soft emission rises from ^1.0 
to ~1.3 keV, and the radius of blackbody slowly increases 
from ~2.5 km to ~2.7 km, the latter values are compat- 
ible to those derived from the subsequent low/hard state 
as reported in Zhang et al. (2009). As the source left the 
quiescent state and started to be active, the cutoff energy 
of the hard component decreased from ~136 keV to ~40 
keV, with the latter being consistent with that derived in 
the subsequent low/hard state as reported in Zhang et al. 
(2009). The confidence contour plots of T and cutoff en- 
ergy are given in Fig. [5] for three INTEGRAL revolutions. 
While some correlation between the two parameters is ev- 
ident in the elliptical shapes of the contours, the different 
fit results are nicely separated, supporting the reality of the 
trend shown in the lower two panels of Fig. |31 



3.2. Type-I X-ray bursts 

111 Chen et al. (2010) we have analyzed 16 bursts that oc- 
c rrred during the 2008 outburst of IGR J17473-2721 and 
found that they could be grouped into two parallel evolu- 
tion tracks regarding their burst duration versus their flux. 
E owever, such parallel tracks did not show up in the work 
by Chenevez et al. (2010), where 57 bursts were found in 
t'le entire outburst, although the overall trend of a positive 
■i;lationship still holds. Instead of taking the 0.1-200 keV 
and by Chenevez et al. (2010) did, we adopt a band well 
wered by RXTE (1.5-30 keV) to minimize the uncertain- 
ties intrinsic to flux extrapolations (see, e.g., Thompson et 
. 2008). 

By using the persistent flux at 1.5-30 keV and the burst 
duration r as derived in Chenevez et al. (2010), we pro- 
duced the r-flux diagram. This is shown in Fig. [6l where 
jse also show linear fits (the parameter pi enclosed in Fig. [5] 
stands for the slope of these linear fits) to the data related 
to bursts occurring in different spectral states. We find that 
the fits worsen when the bursts from HSS and number 8 in 
the classification of Chenevez et al. (2010) are included. 
Recently, Zhang et al. (2010) found while analyzing the 
cooling phase of type-I X-ray bursts in 4U 1636-53 that the 
bursts occurring in HSS may be He-dominated, while they 
could be H-dominated in LHS. Burst number 8 occurred 
less than 10 minutes later after it predecessor, and it most 
likely has the same fuel (Fujimoto et al. 1987). We there- 
fore consider an overall sample for the bursts occurring in 
the LHS that excludes burst number 8. We proceeded to 
check the r-flux relationship in burst groups (the different 
color bars above the PCA lightcurves in the Fig. 1) located 
in the initial rising LHS, the LHS plateau preceding the 
HSS, the whole LHS preceding HSS (i.e., the sum of the 
flrst two groups), the LHS after the HSS, and the overall 
sample (see Fig. ^ . The correlation factors are derived as 
0.90, 0.32, 0.68, 0.75 and 0.35, respectively. A linear fit re- 
sults in a reduced ^ 2.1 for the overall and ~ 1.2-1.8 for 
the rest. These values seem acceptable because the overall 
sample consists of grouped bursts showing different slopes 
in this diagram, for which the is individually better. We 
find that the bursts in the decaying LHS present a linear 
slope of 5.0±1.0, which is substantially different from the 
1.7±0.1 of the overaU and 0.9-1.8 of the rest (see Fig.^. A 
parabolic line slightly improves the fit, with ~62.61 (35 
d.o.f.) for all the bursts except number 8. Fitting all bursts 
except number 8 with a broken line results in a ~53.08 
(34 d.o.f.), and a probability of 0.00296 under an F-test for 
having a marginally-significant turn-over component at 
the lower flux. 

If we enlarge the energy band to 0.1-200 keV, the r- 
flux correlation factors of the four burst groups described 
above are 0.91, 0.29, 0.72, 0.35, and 0.52, respectively. A 
linear slope is derived as 0.7±0.2 for the leading LHS and 
1.9±0.4 for the decaying LHS. Thus, the trend stiff holds 
that the linear slope in the decaying LHS is steeper than 
in the leading LHS, although the persistent flux is esti- 
mated with an energy band (0.1-200 keV) well beyond the 
PCA domain. We also investigate the dependence of ther- 
mal and non-thermal emission, respectively, and obtained 
correlation factors of 0.44 and 0.73 (see Fig. [7]). Because the 
thermal emission is weak in the low/hard state, the error 
bars of the thermal flux turn out to be too large in Fig. 
[71 which prevent us from inferring the r-flux relationship 
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further than a simple overah positive correlation. Moreover 
it is hard to disentangle the disk emission from the thermal 
emission from the neutron star surface. But apart from this 
the diagram of r vs thermal flux as shown in Fig. [7] may be 
used as an standard candle in investigating the dependence 
of the burst duration on the accretion rate. 



4. Discussion 

4.1. Spectral evolution from quiescent state to LHS 

The cutoff energy of the power-law component in X-ray 
spectra is a diagnostic parameter of the electron energy 
distribution in the corona and of the analytical approxi- 
mation of unsaturated Comptonization. The latter is used 
extensively to study the X-ray evolution of BH and NS 
XRB outbursts (e.g., Motta et al. 2009, Chen et al. 2006). 
In general, the spectrum is typical of an LHS with a high 
energy cutoff at ~ 40 keV for the NS binaries and ~ 80 
keV for the BH binaries. However, there are only a few 
sources observed in the X-ray band at the beginning of the 
outburst because either the transition is too quick or the 
source is too faint in this period. This is especially true for 
the NS X-ray binaries because only the low-luminous-type 
atoll sources reach the LHS. For instance, for the LHS of the 
outburst from Aql X-1 in 2005, which has an LHS before 
an HSS analogous to IGR J17473-2721, the cutoff energy 
is ~ 30 keV and no cutoff was detected up to 150 keV in 
the following HSS (Rodriguez et al. 2006). 

In the 2008 outburst of IGR J17473-2721 we have found 
here that the cutoff energy Ecut gradually decreases in 
the state transition from quiescence to LHS state within 
ten days, which is very similar to the spectral evolution 
observed from BH X-ray binaries. Accompanied with the 
spectral evolution shown in our previous work (Zhang et 
al. 2009), evidently sees that i?cut abruptly decreases from 
Ecut = 38.51^gQ^^ keV to '^5 keV when the source moves 
from the LHS plateau to HSS, within three days. For the 
LHS after the HSS, the hard tail of spectra was fitted with 
a bknpow (£:bkn ^^51 keV) or Comptt {Eg ~22 keV). This 
complete evolution of the spectrum follows the sequence of 
quiescent state LHS — ^ HSS — ?> LHS in the outburst of 
IGR J17473--2721. 

From the observations of the transition from quies- 
cent state LHS in IGR J17473-2721 with RXTE and 
INTEGRAL, one sees that the total flux is rising in the 
soft and hard X-ray bands (3-200 keV), the spectra become 
hard, and i^cut decreases monotonously. This resemblance 
in spectral behavior between BH and NS binaries may be 
a clue for understanding the state transitions in both sys- 
tems, which is related to the cooling/heating of the corona 
by the photons produced at the accretion disk or/and the 
NS surface (Esin et al. 1997; Menou et al. 1999). However, 
after the initial LHS, spectral differences between BH and 
NS binaries occur. From our previous work (Zhang et al. 
2009), Ecut of IGR J17473-2721 in HSS was measured to 
be ~5 keV, similar to that of 4U 1608-522 {E^ut --3 keV 
in its HSS, as reported by Chen et al. 2006). These values 
are different from BH systems (where no cutoff is observed 
in HSS). This inconsistency between the two systems may 
be related to the existence of a hard surface or/and the 
magnetic field in NS binaries (Chen et al. 2006, Zhang et 
al 1996). 



4.2. Does most of the accretion material fall to the NS 
surface in the LHS? 

The persistent emission is usually taken as an estimator of 
the accretion rate in XRB; but it is not always easy to find 
out if a significant outflow reduces the final inflow rate onto 
the central compact object. It is also not clear whether the 
inflowing material exhausts all its gravitational potential 
energy before reaching the surface of an NS or the event 
horizon of a BH. Fortunately, under the assumption that 
the burst exhausts the entire accreted material on the NS 
surface, bursts may be used to provide a good estimate for 
the actual inflowing rate; this cannot be done easily for BH 
XRBs. 

The accretion rate estimated from the observed per- 
sistent emission is given by Mp ^ Lp/e, where Lp is 
the persistent luminosity and e is radiative efficiency. In 
NS binaries, the potential energy released per nucleon (of 
mass mp) is GMj^smp/Rj^s ~ 200 MeV (Mns = 1-4 Mg, 
Rns = 10 km), or ^ 20% of its rest mass. We then have 
Mp - 2.4 X 10^"^ g s-^ with Lp ~ 0.11 ^Edd in the LHS 
preceding the HSS of IGR J17473-2721. 

On the other hand, the inflowing rate can be estimated 
from the bursts occurred in the LHS preceding the HSS 
of IGR J17473-2721 by Mb - U/e, where Lb - 1.4 x 
lO"'^ ^Edd is the average luminosity of the bursts in the 
LHS, e ~ Qnuc/938=(1.6-f 4 < X >)/938, Qnuc (in units 
of MeV) is the energy released per nucleon during the burst, 
and X is the hydrogen mass fraction of the fuel. Mb ~ Mp 
requires < X >^ 0.18, a reasonable value and consistent 
with the durations of the observed bursts. We therefore 
suggest that at least for IGR J17473-2721, most of the 
accreted material will fall to the NS surface in the LHS. 

Fig. [7] shows that, similarly to other NS binaries in the 
LHS, the power-law component dominates the persistent 
flux in the two observed LHSs of IGR J17473-2721. This 
means that the majority of the gravitational potential en- 
ergy is released in the corona before reaching the surface of 
the NS; otherwise the kinetic energy of the accreted mate- 
rial would produce luminous thermal (blackbody) radiation 
from the surface of the NS. 

4.3. Spherical or disk corona ? 

The almost pure blackbody spectrum of the observed type-I 
bursts indicates no Comptonization of the thermal photons 
from the surface of the neutron star, i.e., the corona does 
not cover a large fraction of the NS surface. The very differ- 
ent ratios between the observed thermal photons and the 
power-law photons in the quiescent state and in the bursts 
suggest that the soft photons from the surface of NS do not 
dominate the cooling of the corona. An additional evidence 
against spherical corona comes from the previous spectral 
analysis of this outburst (Zhang et al. 2009), where the sys- 
tem was found to be highly inclined and most of the disk 
emission is not supposed to be visible. Additional evidence 
for a disk corona is that a radius ^3 km of the blackbody 
in the persistent flux (the fourth panel of Fig. [5]) is more 
compatible with being a hot spot on the NS surface than 
its inner accretion disk radius. 

On the other hand, a corona on top of the disk has 
a relatively smaller coverage for the soft X-rays from the 
surface of the NS, naturally avoiding all the difficulties dis- 
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Fig. 6. Timescale r vs. the persistent flux in the 1.5-30 keV band for the bursts occurring in LHS. The black points 
represent all of the bursts together (upper left panel). For comparison, each burst group is plotted together with the 
results of a linear fit to this group marked by different colors. The top-left panel shows the fit results for all bursts (with 
larger reduced x^, black) and all bursts except the burst number 8 (in red). The top- right panel shows the fit results of 
the bursts during LHS after the HSS (blue), and the bursts occurring during the rising part of the LHS before the HSS 
(green). The bottom- left panel presents the fit results of the bursts happening at the plateau of the LHS before the HSS 
and the bottom-right panel those for the bursts occurring during the LHS before the HSS (purple), for comparison. The 
fit results for the bursts in the LHS after the HSS (green) are also included. 




flux (10'^ erg cm' s'^) flux (10"' erg crn^ s'^) 

Fig. 7. Plots of the timescale r vs. the thermal (left panel with bolometric blackbody flux) and non-thermal (right panel 
with power-law component in the 1-100 keV band) persistent flux for the bursts occurring in the LHS. 



cussed above. In this case, the non-detection or the low 
level of the disk's thermal X-ray emission would suggest 
that the corona covers the disk almost completely and that 
a saturated inverse Comptonization process produces the 
observed power-law spectrum in the persistent emission. 
We conclude that it is likely that most of the thermal pho- 
tons from the NS surface are not Comptonized, because the 
corona covers the disk, but not the NS. 



4.4. Toward a ms pulsar system? 

By neglecting the influence of a possible outflow on the 
estimation of the accretion rate, assuming a radiative effi- 
ciency of ~ 20% and a typical duty cycle of ^7% (Sivakoff 
et al. 2008), the mass growing-rate of an accreting NS can 
be derived as 



SM 
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Fig. 8. Six-eight percent {la, inner, black), 90% (1.6cr, middle, red), and 99.7%(3cr, outer, green) confidence level contours 
of the high-energy cut-off energy J?cutoff versus the photon index F based on the joint PCA/JEMX/HEXTE/ISGRI 
spectra of rev. 668, rev. 670 and rev. 674 (the last revolution of INTEGRAL in the rising phase) of IGR J17473-2721. 
The best-fit value is indicated by a . 



where Lboi is average NS luminosity during outbursts, e is 
the radiative efficiency, i5 is the duty cycle, and r is the 
lifetime of the NS binary. By taking r ^10^ yr and Lboi ^ 
0.1-lLEddj we find that the NS will accrete a mass of 0.02- 
0.22 Mq from its companion. Millisecond pulsars (MSP) 
with a spin period (Pg <20 ms) seem to have a larger mass 
(Mns=1-57 ± 0.35 Mq) than the mean value of 61 mea- 
sured masses of NSs (1.46 ± 0.30 Mq) (Zhang et al. 2011). 
The estimation above may indicate that a X-ray burster 
such as IGR J17473~2721 could evolve toward an ms-pulsar 
system (Demorest et al. 2010). 



4.5. Type-I X-ray bursts related to the spectral state of the 
outburst 



For the bursts embedded in the outburst, we found that 
for H-dominated burning the CNO cycle constrains the 
speed at which the H fuel is exhausted, which may lead 
to the burst duration correlating with the accretion rate. 
The most recent analysis in Zhang et al. (2010) indicates 
that the bursts at low hard state most likely have H fuel. 
In Fig. ini this correlation is present for the bursts in the 
overall LHS, and the LHS after HSS. The slope of these 
correlations differs at a significance of ~3 a level, with the 
slope for the bursts occurring in the decaying LHS being 
a factor of 3 steeper than that of the overall sample. An 
extrapolation of the slope for the bursts in the decaying 
LHS to the other three groups in Fig. [5] would overshoot 
the lower fluxes. This feature may be related to part of the 
accretion material holding up to form a corona on its way 
of falling onto the surface of neutron star during the de- 
caying part of the low hard state. Therefore, the apparent 
accretion rate on the surface of the neutron star prior to 
each burst would be overpredicted if one simply were to 
take take the persistent emission as estimator. 



4.6. Summary 

With our complementary analysis of INTEGRAL and 
RXTE observations of the 2008 outburst of IGR J 17473- 
2721 we find that 

— The cutoff energy of the spectra decreases from Ecut ^ 
150 keV to ~ 40 keV when the source leaves the qui- 
escent state toward the LHS, which is rarely reported 
in NS binaries. This phenomenon may be related to the 
cooling/heating of the corona by the photons produced 
at the accretion disk or (and or) at the NS surface. 

— The corona seems to be located above the disk, but 
not around the NS, so that emission from type-I X-ray 
bursts can escape from the NS surface without suffering 
Comptonization. 

— The different slopes of the linear relationship hinted in 
the burst duration vs. the persistence flux diagram for 
the LHS before and after the HSS may be caused if part 
of the accretion material is used to form a corona. 
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